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Abstract 

The Selex experiment measured several billion charged hadron tracks with a high precision magnetic momentum 
spectrometer and high precision RICH velocity spectrometer. We have analyzed these data to simultaneously measure 
the masses of all the long lived charged hadrons and anti-hadrons from the n to the Q using the same detector and 
technique. The statistical precision achievable with this data sample is more than adequate for 0.1% mass measure- 
ments. 

We have used these measurements to develop and understand the systematic effects in using a RICH as a precision 
velocity spectrometer with the goal of measuring 10 masses with precision ranging from 100 KeV for the lightest 
to 1000 KeV for the heaviest. This requires controlling the radius measurement of RICH rings to the ~ 10~ 4 level. 
Progress in the mass measurements and the required RICH analysis techniques developed are discussed. 



The Selex RICH was orginally conceived as a parti- 
cle identifier for a fixed target multi-particle spectrom- 
eter Once we saw real data from this detector we 
learned that this technique had serious potential as a 
precision velocity spectrometer J2> S H S] ■ Plots like 
Figure Q] ||2|] clearly demonstrated that precision mass 
measurements of many of the hadrons, at the same time 
with the same detectors, was possible. This paper is an 
exploration of the systematic limits of this technique. 

Recently the MIPP experiment proposed further data 
taking to resolve the ~ 100 ppm discrepancy in the latest 
charged kaon x-ray mass measurements [6]. They now 
have the Selex RICH [Q] (with C0 2 instead of Ne as a 
radiator) but with insufficient data to reach the 100 ppm 
level of statistical precision. Selex has large amounts 
of data already recorded, and analysed, in a well under- 
stood apparatus. Selex events are multi-hadronic in- 
teractions, not single tracks. Nonetheless we thought it 
would be good to see how far we could take the idea 
of making precision particle mass measurements using 
a precision momentum (magnetic) and velocity (RICH) 
spectrometers. We have more than enough data to make 
mass measurements with statistical precsion better than 
the < 1000ppm= 0.1% level, even with very tight cuts 
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to select clean individual track measurements. 

The goals of this study are to identify the important 
systematic uncertainties of RICHs used a precision ve- 
locity spectrometers. Only studies with real data can 
fully illuminate the the systematic limitations of such 
precision spectrometry. It worth noting that m - p/(J3~y), 
where p is the momentum, and f3y the relativisitic veloc- 
ity, so that systematics in the magnetic momentum spec- 
trometry are equally important: both must be in control 
to achieve a precision mass measurement to the < 100 
ppm level. 

Selex was a Fermilab fixed target experiment de- 
signed to study the production and decay of charmed 
baryons. It took data in 1996 - 7 in a 600 GeV/c E" 
beam with excellent vertex and momentum spectrome- 
ters. The Selex RICH, one of the first large multi-pixel 
PMT RICHes, provided ~ 1% velocity resolution for 
all particles above threshold (J3y > 86). This is well 
matched to the 0.5 - 1 % momentum resolution for tracks 
which reach the RICH. 

The same detectors and analysis provide common 
systematics for all particles species. Figure [1] displays 
18 particles species from a 12. 5M track sample with 
RICH Ring Radius plotted as a function of measured 
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track momentum. We have no sensitivity to the elec- 
tron and muon masses and we can't resolve the S + from 
the anti-2~T. We can measure the masses of the other 10 
particles. 



mass values are plotted. The statistical mass uncer- 
tainties are (n = 70 KeV/c 2 , K = 160 KeV/c 2 and 
p = 350 KeV/c 2 ). Systematic uncertainties clearly 
dominate these mass measurements. 




100 200 300 

P [GeV/c] 



Figure 1: RICH Ring radius vs measured momentum for 12. 5M 
tracks. 



The usual small angle and ultra-relativistic approx- 
imations [& c = ^26 — 1/y 2 ] to the Cherenkov equa- 
tion [cos(& c ) = l/0n 6 = n - 1 = 67.xl0" 6 ] is good 
to 0.05%. The measured ring radius is related through 
the spherical mirror focal length (F) to the momen- 
tum, mass and maximum Cherenkov angle by R(p) - 
FyJ(&'? ax ) 2 -(m/p) 2 . 

Exploiting this relationship to measure mass requires 
three calibrations; the momentum scale, determined by 
reconstructing K 1 ^ decays, the mirror focal length, which 
varies from 989 - 992cm across the 16 mirrors in the 
Selex RICH, and the maximum Cherenkov angle 0" ,flA . 

Rings are fit for the radius and center coordinates |0] 
from a PMT list generated around the measured track 
angles in the RICH cut around the predicted radius for 
a given mass hypothesis. At least 5 hits are required. 
As a clean place to start the Selex data used are low 
statistics, low rate, low multiplicity data with one of the 
two spectrometer magnets off. A mass spectrum gen- 
erated after an initial calibration for tracks with ring 
radii in the interval 5 - 8cm (R max = 11 Acrri) is shown 
in Figure |2|top). Gaussian mass fits in 2mm R bins 
give poor^ 2 mass averages shown in Figure ^bottom, 
open black points). The difference with the PDG 
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Figure 2: Mass from measured ring radius and momentum [top]. Fits 
to masses in 2mm radius bins; open (before), and solid (after) accep- 
tance corrections [bottom] . 



The momentum resolution is shown in Figure [3] For 
the data set analysed here the momentum cutoff for a 
track to reach the RICH is 12 GeV/c. For the Selex 
charm data the resolution is much better but the cutoff 
is 20 GeV/c, putting the pions out of reach but greatly 
improving the momentum resolution for the hyperons. 
The momentum scale was calibrated using — » jt + jt~ 
decays and the mass from the PDG [9]. With this 
calibration the masses reconstructed with charged parti- 
cle decays of 12 other hadrons from the <p through the 
Q.® are correctly reproduced. 

We've identified several systematic effects. The first 
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is just geometry; as shown in Figure [4| top), the inter- 
section of a circular ring with a circular tube isn't sym- 
metric. The acceptance correction (proportional the the 
arc length shown) is biased towards radii larger than the 
center of ring to center of tube distance. This accep- 
tance is a just geometric. It is plotted as a function of 
z = R' — R in Figure |4f bottom) for different ring radii. 

We've build a second ring radius fitter based on max- 
imizing the joint likelihood for all tubes on a ring 
where the likelihood for each tube is from the accep- 
tance curves shown. The average shift in ring radius is 
~ 0.1mm or 0.1 - 0.2% for rings in the radius interval 
(5-8 cm) we use to fit the mass. 

This likelihood fitter has an interesting application in 
pattern recognition. The acceptance plotted goes ex- 
actly to zero at the radius difference where the ring no 
longer crosses the hit tube. To deal with these cases 
we assign a noise probability of 0.002 for the minimum 
likelihood of any tube in the ring fit. The likelihood 
as a function of fit radius can have one or more max- 
ima. Consider the case of two tubes along a ring radius: 
only one can fit the ring, the other being "noise" in this 
model. This likelihood has two maxima (if all the other 
tubes in the fit are well behaved). Rejecting these am- 
biguous fits is an objective way to remove non-Gaussian 
tails due to pattern recognition mistakes from the ring 
radius resolution function. 
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Figure 3: Momentum resolution for the data set analysed here (solid) 
and for the Selex charm data (dashed). The momentum region corre- 
sponding to 5 < R(an) < 8 for each particle species is also shown. 



Repeating the mass analysis after applying cuts based 
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Figure 4: The intersection of a ring with a tube isn't symmetric. The 
acceptance is biased to R beyond the tube center [top]. The acceptance 
function as a function of ring radius [bottom]. 



on the criteria outlined above yield the mas values 
shown in Figure [2] The x 1 for the mass averages as 
a function of radius are improved but systematics in 
the mass as a function of ring radius remain for the pi- 
ons and kaons when the calibration constants are deter- 
mined using the protons. 

We have made some progress. The proton and anti- 
proton mass now agree with each other. We used the 
protons for calibration so the proton mass has to agree 
with the PDG value. The kaon masses are close to each 
other and to the PDG values but still several 100 KeV 
and several <x off. The pion masses are still quite far 
from the PDG values for reasons which are unclear. We 
can and will apply these methods to the charged hyper- 
ons when the systematics are better understood. 

This study is beginning to illuminate some of the sys- 
tematics of the resolution of RICHs as precision veloc- 
ity spectrometers. More work will be required to reach 
and understand the resolution limits of this technique. 

We are indebted to Selex for the RICH data, our 
home institutions, Consejo Nacional de Ciencia y Tec- 
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nologfa (CONACyT), Fondo de Apoyo a la Investi- 
gacion (UASLP), and the U.S. Department of Energy 
(contract DE-AC02-76CHO3000), for support. 
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